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ABSTRACT 

Aims. We compare the performance of several dust models in reproducing the dust spectral energy distrihution (SED) per unit extinc¬ 
tion in the diffuse interstellar medium (ISM). We use our results to constrain the variability of the optical properties of big grains in 
the diffuse ISM, as published by the Planck collaboration. 

Methods. We use two different techniques to compare the predictions of dust models to data from the Planck HFl, IRAS and SDSS 
surveys. First, we fit the far-infrared emission spectrum to recover the dust extinction and the intensity of the interstellar radiation field 
(ISRF). Second, we infer the ISRF intensity from the total power emitted by dust per unit extinction, and then predict the emission 
spectrum. In both cases, we test the ability of the models to reproduce dust emission and extinction at the same time. 

Results. We identify two issues. Not all models can reproduce the average dust emission per unit extinction: there are differences of 
up to a factor ~ 2 between models, and the best accord between model and observation is obtained with the more emissive grains 
derived from recent laboratory data on silicates and amorphous carbons. All models fail to reproduce the variations in the emission 
per unit extinction if the only variable parameter is the ISRF intensity: this confirms that the optical properties of dust are indeed 
variable in the diffuse ISM. 

Conclusions. Diffuse ISM observations are consistent with a scenario where both ISRF intensity and dust optical properties vary. The 
ratio of the far-infrared opacity to the V band extinction cross-section presents variations of the order of 20% (40 - 50% in extreme 
cases), while ISRF intensity varies by ~ 30% (~ 60% in extreme cases). This must be accounted for in future modelling. 


1. Introduction 


Interstellar dust is an important component of the interstellar 
medium (ISM): it is studied for its role in the physics and chem¬ 
istry of the ISM (D’Hendecourt et al. 1985; Mathis 1990; Bakes 
& Tielens 1994; Le Bourlot et al. 2012), for its role as a tracer of 
gas column density (Boulanger et al. 1996; Rowles & Froebrich 
2009; Planck Collaboration XI 2014) and of magnetic field 
structure (Chapman et al. 2011; Poidevin et al. 2013; Berdyugin 
et al. 2014; Planck Collaboration Intermediate Results XIX 
2015), and for its effect as a foreground in studies of the cosmic 
microwave background (Planck Collaboration Int. XXII 2015; 
Planck Collaboration Int. XXX 2014). Many dust models have 
been made to reproduce the main dust observables, which in¬ 
clude extinction curve and albedo, spectral energy distribution 
(SED) from the near-infrared to the microwave continuum, and 
elemental abundance constraints (Zubko et al. 2004; Draine & 
Li 2007; Compiegne et al. 2011; Jones et al. 2013). 

Our understanding of interstellar dust, however, is still in¬ 
complete. One of the issues faced by the models is that dust is not 
the same everywhere in the ISM. As revealed by numerous re¬ 
cent observations (Pagani et al. 2010; Kohler et al. 2011; Planck 
Collaboration et al. 201 lb, a; Kohler et al. 2012; Martin et al. 
2012; Paradis et al. 2012; Ysard et al. 2013), dust properties vary 
from the diffuse ISM to molecular clouds. Variations are espe¬ 
cially notable in the far-infrared opacity, i.e. the far-infrared op¬ 
tical depth per unit dust column density tINh. The dust optical 
depth is the function that, in the optically thin limit, modu¬ 


lates the SED: h - tx ■ Bx(T), where B,i(T) is the black-body 
emission at a temperature T * . 

Variations in dust far-infrared opacity can now be traced in 
the diffuse ISM itself thanks to the survey of the Planck satel¬ 
lite (Tauber et al. 2010). The far-infrared opacity shows signifi¬ 
cant variations in the diffuse ISM, although those variations are 
smaller than the variations observed between the diffuse ISM 
and molecular clouds (Bot et al. 2009; Planck Collaboration Int. 
XVII2014; Planck Collaboration XI 2014). 

Planck Collaboration Int. XXIX (2014) examines the dust 
emission obtained from Planck , IRAS (Neugebauer et al. 1984) 
and WISE (Wright et al. 2010) data, fitting it over the whole sky 
with the Draine &. Li (2007) dust model, which uses constant 
optical properties ^ . In this paper, our analysis concentrates on 
more than 200,000 lines of sight towards quasi-stellar objects 
(QSOs) observed in the Sloan Digital Sky Survey (SDSS, York 
et al. 2000) and for which Ay is available ^ : the Ay from the 
QSOs is correlated with the value A^ deduced from the Planck 

* Quantities that depend on dust optical properties, such as optical 
depth or opacity, are wavelength-dependent. When, in this article, we 
mention these quantities without specifying a wavelength, the quantities 
should be interpreted as an “effective” value integrated over the whole 
wavelength range of interest. For instance, since we are interested in 
the far-infrared emission of grains, we often mention the “far-infrared 
opacity”: this is the integral of the opacity over the far-infrared and 
submillimeter wavelength range, weighted by the emission spectrum. 

^ The optical properties for the materials used in Draine & Li (2007) 
were modified as explained in Aniano et al. (2012). 

^ The SDSS survey actually provides the reddening E(B - V), from 
which Ay is computed assuming Av/E(B - V) = 3.1. Throughout the 
present paper, as in Planck Collaboration Int. XXIX (2014), when- 
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and IRAS fit. Planck Collaboration Int. XXIX (2014) find that 
Galactic dust SEDs can be well fitted by the Draine & Li model 
(Draine & Li 2007) when one uses as free parameters the dust 
column density - or, equivalently, the V band extinction Ay - 
and the intensity Umin of the interstellar radiation field (ISRL) 
in the diffuse ISM. However, even in the diffuse ISM, this model 
cannot explain the variations in the SED per unit extinction if the 
ISRL intensity is the only quantity that varies : the Draine & Li 
(2007) model can reproduce the dust SED, but it systematically 
overestimates the observed Ay by a factor that depends on f/njin 
and has a mean value of 1.9 in the diffuse ISM. 

We aim to extend the analysis by Planck Collaboration 
Int. XXIX (2014) to the Compiegne et al. (2011) and Jones 
et al. (2013) dust models and to quantify variations of their far- 
infrared emission properties in the diffuse ISM. Our analysis will 
focus on the big grains in thermal equilibrium (radius ~ 100 nm). 
The paper is organized as follows. In Sect. 2 we introduce the 
observational data and the dust models used to fit these data. In 
Sect. 3 we present two different methodologies for constraining 
dust parameters using observations, we show how their results 
differ and we estimate the variation of dust optical properties in 
the diffuse ISM. In Sect. 4 we discuss the results, identifying 
a possible shortcoming of most current dust models. Linally, in 
Sect. 5 we explore the prospects for future dust models. 

2. Data and models 

2.1. Data 

Planck Collaboration Int. XXIX (2014) provides a family of 20 
Ay-normalised SEDs, Ti/Ay, to be used as constraints for fu¬ 
ture dust models. The quantity Ay is measured on QSO lines- 
of-sight: each SED is the average over ~ 10,000 lines of sight 
with similar U^m- The SEDs have t/min in a range of 0.42 to 
0.98, with an average of 0.66 and a standard deviation of 0.14. 
The QSO lines of sight sample regions of low Ay with a median 
value of ~ 0.10, however, the set of SEDs describe the data over 
a much larger fraction of the sky (> 70%) with Ay up to 1 The 
wavelengths coverage is from 60 /rm to 2.1 mm (between 5000 
and 143 GHz). 

We use a selection from the Planck Collaboration Int. XXIX 
(2014) SEDs. Dust models agree in attributing the V band ex¬ 
tinction mainly to big grains, which, being in thermal equilib¬ 
rium, also dominate emission at T > 100 pm. To compare the 
same grain type in extinction and in emission, therefore, we 
use the 100 pm (3000 GHz) IRAS band and the 350, 550, and 
850 pm (857, 545 and 353 GHz) Planck bands. We exclude 
the 60 pm band, where emission is dominated by stochastically- 
heated smaller grains, because a single band is not sufficient to 
test small grains properties; a measure of the extinction in the 
UV is needed to determine the total mass of small particles, and 
a full mid- and near-infrared emission spectrum is needed to con¬ 
strain the small particle temperature and size distribution. 

Eig. 1 shows a subset of the data; the median SED, with 
Umm - 0.66; the two extreme SEDs, with t/min = 0.98 and 
Urain - 0.42; the SEDs closest to being 1 cr above and below 
the median, with t/min = 0.80 and 0.52 respectively. The figure 
shows that the Ay-normalised SEDs are very similar over the 

ever the observed Ay is mentioned, the measure referred to is actually 
3.1 • E(B-V). 

The advantage of normalising dust emission to Ay instead of the HI 
column density is that the data contain information on dust properties 
alone, without the uncertainty introduced by the dust-to-gas mass ratio. 

^ G. Aniano, private communication. 



X (/Am) 

Fig. 1. Subset from the 20 Ay-normalised dust SEDs in Planck 
Collaboration Int. XXIX (2014). Each Umm, i e. each tempera¬ 
ture, is shown in a different colour and line style. Eor clarity of 
illustration, only five of the 20 SEDs are shown; the warmest and 
coldest overall SEDs, the warmest and coldest of the SEDs that 
are within a standard deviation from the average Umm, and the 
median SED. Symbols indicate the central band wavelengths; 
error bars indicate the dispersion for each band flux, while in¬ 
strumental noise is negligible. 

Planck wavelength range, and most of the variations are in the 
IRAS 100 pm band. 

The SEDs have both statistical and systematic uncertain¬ 
ties. What is important in statistical uncertainties is their rela¬ 
tive value from one band to another, which defines the statis¬ 
tical weight to be applied to each band in the fitting routine. 
The statistical errors in each band of /,i/Ay are negligible, since 
each SED is an average over more than 10,000 observations. 
The systematic uncertainties in the dust SED come mainly from 
the photometric calibration; the calibration uncertainty is impor¬ 
tant at 100 pm (Icr error bar = 13.6%) and at 350 and 550 pm 
(Icr = 10%); we decided to neglect it at longer wavelengths 
where it is much smaller (^ 1%) (Miville-Deschenes & Lagache 
2005; Planck Collaboration VIII 2014). We estimated the effect 
of this uncertainty via Monte Carlo simulations. Specifically, 
for each observed SED, we realized 1000 simulations of ran¬ 
dom Gaussian-distributed errors for the 100, 350, and 550 pm 
bands, with the photometric error at 350 and 550 pm set equal 
since both channels are calibrated on planets (the uncertainties 
are dominated by the planet model used, not by noise). We per¬ 
formed a fit on all simulations, obtaining 1000 values of each 
fit parameter for each SED; the mean of these was taken as the 
fiducial value of the fit parameter and the standard deviation as 
the uncertainty. The systematic uncertainty in Ay is estimated by 
Planck Collaboration Int. XXIX (2014) to be ~ 15%. 

2.2. Dust models 

We use three dust models with different characteristics; Draine 
& Li (2007), with the modifications explained in Aniano et al. 
(2012), Compiegne et al. (201 1), and Jones et al. (2013). 

Draine & Li (2007), hereafter DL07, is one of the most pop¬ 
ular dust models, often used to estimate dust and gas masses in 
Galactic and extragalactic environments by fitting dust emission 
(Draine et al. 2007; Aniano et al. 2012). It comprises a size dis¬ 
tribution of big silicate grains with optical properties adapted 
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from observations (Draine & Lee 1984; Li & Draine 2001) and 
carbonaceous grains whose optical properties are assumed to 
vary continuously from those of graphite for large sizes (» 10^ 
C atoms) to those of polycyclic aromatic hydrocarbons (PAHs) 
for small sizes 10^ C atoms). 

Compiegne et al. (2011), hereafter Cll, uses silicate grains 
with the same properties as DL07, but carbonaceous grains in 
Cll are divided into two populations: one of PAHs, and one 
population of amorphous carbon grains, covering a size range 
from nm-sized to big grains, with optical properties derived from 
Zubko et al. (1996). 

Jones et al. (2013), hereafter J13, is a core-mantle model 
with optical properties that are derived from laboratory mea¬ 
surements. The J13 model features big grains of amorphous 
forsterite-like silicates, with iron inclusions and amorphous 
aromatic carbon mantles; big amorphous carbon grains with 
aliphatic cores and aromatic mantles; and small aromatic grains 
(down to sub-nm sizes), which take the place of PAHs. A unique 
feature of J13 is that it allows for a certain flexibility in its opti¬ 
cal properties, which vary according to the band gap of carbona¬ 
ceous materials and the size of small grains, as well as possible 
variations in the mantle thickness. We use the standard model 
aromatic carbon with a band gap of 0.1 eV. 

The dust extinction and emission are computed by interpo¬ 
lating an existing library in the case of DL07 (Aniano et al. 
2012) and using the DustEM tool ^ in the case of Cll and J13 
(Compiegne et al. 2011). The observational SEDs were fitted us¬ 
ing a^^-minimising procedure. We assumed an interstellar radi¬ 
ation field (ISRE) of the form M,i = where is 

the mean ISRE in the solar neighbourhood estimated by Mathis 
et al. (1983) and Go is a dimensionless scale factor. The factor 
Go has the same meaning in the model as the Umia introduced in 
Sect. 1 , but we chose to use a different name to more easily dis¬ 
tinguish our results and those of Planck Collaboration Int. XXIX 
(2014). Eor the rest of the paper. Go is our estimate of the ISRE 
intensity, while Gmin is only used to indicate the SEDs and is 
functionally synonymous with dust temperature. 

3. Methodology and results 

We analysed the observations using two different methods. In 
the first method, we fit the SED to obtain the model dust param¬ 
eters: the radiance, i.e. the integrated SED intensity, the ISRE 
intensity Go, and the dust column density, as represented by Ay. 
In the second method, we estimate the ISRE intensity from the 
observed Ay and dust radiance, and then use it to predict the 
model dependent SED. 

3.1. SED fitting 

The ISRE intensity Gq‘ is obtained from the fit of the SEDs. The 
same fit, since the SEDs are normalised by the QSO-derived Ay, 
returns the ratio of the model extinction to the observed extinc¬ 
tion, A^'/Ay. The fit results are shown in Eig. 2; each symbol 
represents a pair of values (Gq‘, A®7^v) for a different model 
and SED. A perfect model would find a value of Ay /Ay - 1 
(thick horizontal grey line). The thinner horizontal grey lines 
show our 15% systematic uncertainty in Ay (see Section 2.1). 
This uncertainty is identical for all SEDs, and so it affects the 
significance of the average Ay /Ay of the models, but does not 
affect either the dependence of Ay /Ay on Go or the differences 

® http://www.ias.u-psud.fr/DUSTEM/ 
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Fig. 2. Plot of the fitted parameters Gjj' and Ay/Ay, for the three 
models. The colour scale goes from blue (cold) to red (warm). 
The filled symbols correspond to the SEDs shown in fig. 1. All 
points have similar relative uncertainty; the error bars for the 
lowest Gmin are shown. The uncertainties are systematic, so they 
affect each point in the same way; errors on the two axes are 
strongly anti-correlated (see Sect. 3.1). The horizontal grey lines 
show A^/Ay - 1 and the 15% uncertainty in the Ay normalisa¬ 
tion (Sect. 2.1). The slanted grey lines are the least-square fits for 
each model; they were added to more easily tell apart the three 
models. 

between models. The error bars on (Gq^ A^/Ay) are obtained 
from Monte-Carlo simulations, as explained in Sect. 2.1. They 
are strongly anti-correlated, with an average Pearson correlation 
coefficient of -0.77 (DL07) to -0.84 (Cll and J13). Since these 
uncertainties are systematic, the plots may be shifted vertically 
or horizontally, but their shape remains essentially the same. 

The J13 model fits the data well, with an average A^ which 
coincides with the expected value. The Cll model overestimates 
Ay by ~ 25%. The DL07 model shows the largest discrep¬ 
ancy, with an Ay that is overestimated by a factor ~ 1.8, as 
already pointed out by Planck Collaboration Int. XXIX (2014). 
All models, however, show a negative correlation between Go 
and Ay /Ay whereas Ay /Ay should be unity ^ . Eor comparison, 
we also fitted the data with a modified black-body. We obtain 
similar results, which are presented in Appendix A. 

3.2. Recovery of Gq from dust radiance 

The SED fitting yields biased results if the optical proper¬ 
ties of the model are not optimal, as was shown by Planck 
Collaboration Int. XXIX (2014) for the DL07 model. It would 
be useful to find an alternative way to compare models and ob¬ 
servations. Therefore, we decided, in addition to fitting the dust 
SEDs, to recover Gq using a procedure based on the dust radi¬ 
ance per unit extinction. 

The radiance, 7?, is the total power emitted by dust grains in 
thermal equilibrium (Planck Collaboration XI 2014). Through 
the conservation of energy, 7? corresponds to the power absorbed 
by grains in thermal equilibrium and can be considered a “heat- 

^ While the trends for Ay and Go in Fig. 2 resemble the well-known 
T-t anti-correlation caused by noise (Shetty et al. 2009a, b), this is not 
the case here. Since each SED is an average over 10,000 spectra, the 
statistical uncertainty for each SED is negligible. 
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Fig. 3. Estimation of Go from the interpolation of 'RlAy. The 
horizontal lines are the value of fi/Ay for the observed SEDs, the 
curves from the origin are the The abscissae 

of the intersection points are the estimates for Go for that SED 
and model. The same 5 SEDs as Eig. 1 are shown, for clarity of 
illustration, using the same colour and line style scheme. 


ing power” for the dust, which is only determined by its absorp¬ 
tion properties, independent of its far-infrared opacity. 

The thermal emission of big grains can be reasonably well 
fitted with a modified black-body, - B^{T) ■ tq ■ {Aj with 
T the observational dust temperature, jS its spectral index and tq 
its optical depth at the reference wavelength To (here we use 850 
yum, equivalent to a frequency vq = 353 GHz). The integration 
of a modified black-body yields (Planck Collaboration XI 2014) 


TT i/rvoj r(4)^(4) 


( 1 ) 


The functions T and ^ are the Gamma and Riemann zeta func¬ 
tion, respectively, and cr^ the Stefan-Boltzmann constant. 

Eor a given dust model one can compute its Ay-normalised 
fluxes in the IRAS and Planck bands, fit a modified black-body 
to the normalised SED (for T > 100 jum), and substitute the re¬ 
sulting parameters in Eq. 1. The result, model^ ^ Pnear 

function of Go, and can be inverted to find Go from the radi¬ 
ance of an Ay-normalised SED. Eig. 3 shows 'R™°‘*‘^'(Go)/A™°‘*‘^' 
for the three models. As we mentioned, "R only depends on ab¬ 
sorption, and so the estimation of ISRE intensity given by this 
method, G^, is not affected by variations in dust opacity. Since 
the method of obtaining Gq assumes that the dust optical prop¬ 
erties are fixed in extinction, while the method of obtaining G|j‘ 
assumes that optical properties are fixed in both extinction and 
emission, we expect Gq to be much less biased than G^f This 
conclusion is supported by the fact that Go estimates become 
much more model dependent when they include optical prop¬ 
erties in extinction: Gg varies by ~ 15% between models, as 
opposed to the factor ~ 2 in Gg'. 

The values for Gg, shown in Eig. 3, vary by a factor 1.6 be¬ 
tween the two extreme SEDs and by 1.3 if we only consider the 
SEDs within Icr of the average (which would represent the “typ¬ 
ical” Go variation rather than the full range). In comparison, in 
Eig. 2, where the dust optical properties in the model are fixed 
and the ISRE is the only source of difference between SEDs, 
Gg‘ varies by a factor 2.3 between extremes and 1.6 within the 
“+lcr” range. 


After obtaining Gg, we use it as a fixed parameter in the 
computation of the SED of the three dust models. The results are 
shown in Eig. 4. To highlight spectral shape variations, the SEDs 
have been normalised by the median data SED and then further 
divided by their RjAy. The J13 model reproduces the median 
SED almost perfectly. Predictions based on the Cl 1 model sys¬ 
tematically present an excess emission at short wavelengths and 
low emission at long wavelengths. This situation is even more 
drastic in the DL07 model. Eurthermore, for all models, the nor¬ 
malised SEDs in Eig. 4 spread over a smaller range than the data: 
the warm (high Gmin) SEDs and the cold (low Gmin) ones are 
closer in temperature for the model than they are for the obser¬ 
vations. 

The differences shown in Eig. 4 between the J13, Cll and 
DL07 dust models can be attributed to their optical properties 
in far-infrared emission, as opposed to their optical properties in 
extinction We find that the J13 model (Jones et al. 2013) is 
a better model for the diffuse ISM emission than that presented 
in Compiegne et al. (2011) and an even better model than the 
Draine & Li (2007) model. The J13 model uses optical proper¬ 
ties for silicates and amorphous carbon based on lab measure¬ 
ments. The Compiegne et al. (2011) and Draine & Li (2007) 
models use silicates with far-infrared optical properties extrap¬ 
olated from mid-infrared astronomical observations (Draine & 
Lee 1984) and empirically adjusted (Li & Draine 2001) to repro¬ 
duce the LIRAS spectra of Einkbeiner et al. (1999) these opti¬ 
cal properties pertain to silicates that are contaminated by other 
materials, including possibly a carbonaceous dust component. 
Eurthermore, Draine & Li (2007) is unique among the models in 
that its carbon dust is assumed to be graphite. Apparently, “as¬ 
tronomical silicate” and graphite are not emissive enough in the 
far-infrared and submillimeter. 

3.3. Variation in dust opacity 

The models shown in Eig. 4 assume constant optical proper¬ 
ties and are tuned to reproduce the correct R/Ay for each SED, 
meaning that spectral shape variations in the modelled SEDs are 
solely driven by the ISRE. It is clear that the full range of spectral 
shapes exhibited by the observations (Eig. 4, top left) cannot be 
reproduced by ISRE variations alone. As mentioned in Section 
3.2, if the dust optical properties are fixed the variations of Gg 
are only about half of what is needed to reproduce the data. We 
need dust with variable optical properties to reproduce the ob¬ 
servations. 

We start by analysing a simplified scenario in which the dust 
optical properties in extinction are fixed, i.e. Ay is a proxy for the 
dust column density, and the SED shapes are explained by the 
variation in far-infrared opacity. Our aim is to quantify this vari¬ 
ation, at least to a first approximation. At the end of this section, 
we see what changes when the assumption of constant extinction 
is dropped. 

Since the SEDs in our work are normalised by Ay, any deter¬ 
mination of the far-infrared optical depth from the data actually 
yields TpiR/Ay. In addition, since Ay is a proxy for the dust col¬ 
umn density, its variations are identical to those of the opacity. 
Eig. 2 shows that Ay/Ay increases by 40 - 50% between the 
warmest and coldest SEDs, independent of the model. This in- 


® The models have very similar optical properties in extinction, see 
Section 3.2. 

^ The high-Galactic latitude Finkheiner et al. (1999) spectra, extrap¬ 
olated from IRAS, predicted a different SED than that subsequently ob¬ 
served by Planck (Planck Collaboration Int. XXIX 2014). 
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Fig. 4. Data SEDs (top left) and dust model I^I^R, normalised by the median data SED (t/min = 0.66) to highlight the spectral shape 
variations. The model SEDs are computed using Gq as ISRE intensity (Sect. 3.2). The same five SEDs as Eig. 1 are shown, for 
clarity of illustration, using the same colour and line style scheme. The observations are plotted in grey behind the models to aid 
comparison 


dicates that in the diffuse ISM, tfir/Av (and therefore opacity) 
varies by a factor 40-50% between extremes. When considering 
just the SEDs within Icr of the average, we obtain an estimate of 
the “typical” variation; ~ 20%. 

We can compare this estimate to another, obtained with a dif¬ 
ferent technique. We repeat the analysis of Section 3.2 on mod¬ 
ified C11 and J13 models in which the far-infrared opacity of 
silicates and carbon has been multiplied by the same factor, in¬ 
dependent of grain size, as shown in Eig. 5. This modification 
is a purely phenomenological artifice, not based on any phys¬ 
ical properties of the materials. We only use this modification 
to estimate the variations in far-infrared opacity that a physical 
model should attempt to reproduce. We find that Cl 1 would re¬ 
produce all the SEDs within Icr of the average if we allowed its 
far-infrared opacity to vary between 1.1 and 1.3 times its stan¬ 
dard value (between 1.0 and 1.4 to reproduce the full range of 
observations). The results are similar for the J13 model with 
scaling factors 0.9 and 1.1 for the Icr range of SEDs (0.8 and 
1.2 for the full range). Thus, we estimate the variations of far- 
infrared opacity to be ~ 20% (40 - 50%) for the typical (full) 
range of diffuse ISM. 

So far, we have assumed that optical properties in extinc¬ 
tion remain the same while the opacity varies. If we drop the 


assumption the conclusion must be modified. When the optical 
properties in extinction vary. Ay is no longer a proxy for the 
dust column density: our previous conclusions on the variations 
of Tfir/Av are still valid, but they no longer apply to opacity 
either. Since Ay is proportional to the extinction cross-section 
in the V band, this means that the variation range found above 
applies to the ratio of the far-infrared opacity to the extinction 
cross-section in the V band. 


4. Discussion 

While some models are closer to the mean SED, no model is able 
to reproduce the full extent of observed variations in the SEDs. 
We recall that in the three models the optical properties of dust 
are assumed constant. Our results show that the variation of Go 
alone is not sufficient to reproduce the entire span of spectral 
shapes presented by the observations. Models that endeavour to 
explain the emission of the diffuse ISM need to include physical 
processes that change the far-infrared optical depth per unit Ay, 
or the ratio of the far-infrared opacity to the V band extinction 
cross-section. 

Variations in size distribution can lead to changes in Ay 
while having little effect on the far-infrared opacity. However, 
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Fig. 5. Effect of variable far-infrared opacity on Z,)/?? for the Cll (left) and J13 (right) models. The same five SEDs as Eig. 1 
are shown, except the median SED, using the same colour and line style scheme. The normalisation is the same as Eig. 4. The 
observations are plotted in grey behind the models; the standard dust models are the solid lines with filled polygonal symbols and 
the modified dust models are the dashed lines with empty symbols. Numbers show the multiplicative factor of far-infrared opacity. 


the resulting variations in 4;/Ay are small given the constraints 
imposed on Go by ’RjAy (See Appendix B). A promising route is 
to consider variations in the grain composition. One process that 
could alter the grain composition is the accretion of carbon, from 
the gas phase or from small particles, on big grains. It is readily 
possible to consider this process within the J13 model because 
it is tuned to exactly this kind of scenario; its big grains, both 
silicate and carbonaceous, include an aromatic carbon mantle, 
the thickness of which can be modified. In the same model, it is 
also possible to explore the effect of having iron inclusions par¬ 
tially composed of iron sulphide (Kohler et al. 2014). Ysard et al. 
(2015) shows that the two scenarios, together with variations in 
the grain size distribution and in the relative abundance of sili¬ 
cate and carbon grains, can explain most of the observed SED 
variations. Grain shape and structure also play a role in opac¬ 
ity. It is known, for instance, that including porosity in grains 
increases the far-infrared opacity, while leaving the properties 
in extinction relatively constant (Ossenkopf 1993; Stepnik et al. 
2003; Kohler et al. 2012). Porous grains are usually considered 
typical of the dense ISM rather than the diffuse ISM, but it has 
been suggested (e.g. Martin et al. 2012) that they could migrate 
from the dense environment to the diffuse ISM. Ear-infrared 
opacity also increases when one employs non-spherical grains 
(e.g. Siebenmorgen et al. 2014), the existence of which is at¬ 
tested by interstellar polarisation. Determining whether a rela¬ 
tion between dust evolution and polarisation exists may be an 
interesting follow-up analysis of the data. 

5. Conclusion and perspectives 

We confront three dust models (Draine & Li 2007; Compiegne 
et al. 201 1; Jones et al. 2013) to a family of Ay-normalised SEDs 
obtained from a combination of the Planck , IRAS, and SDSS 
surveys probing the Galactic diffuse ISM (Planck Collaboration 
Int. XXIX 2014) . 

The widespread method of fitting the SED to obtain the ra¬ 
diation field intensity Go and the dust column density (from Ay) 
gives results that are very model dependent. The average A^ of 
Jones et al. (2013) is close to the observed value, while that of 
Compiegne et al. (201 l)is about 25% too high and that of Draine 
& Li (2007) too high by a factor ~ 1.8. Eurthermore, all mod¬ 


els present a non-physical mass-temperature anti-correlation and 
yield biased Go estimates. 

We develop an alternative Go estimate that uses the radiance 
per unit extinction, 'RfAy. This estimate is robust with respect 
to variations in the far-infrared optical properties and gives less 
biased results. We use these values to compute the SEDs for the 
models, finding that the variations in observed SEDs cannot be 
reproduced without varying dust optical properties. We estimate 
that, to reproduce the typical variations observed, we need to 
vary the ratio of the far-infrared opacity to the V band extinction 
cross-section by ~ 20% in the typical diffuse ISM (up to 40 - 
50% to explain the full range of observations); at the same time, 
variations of ~ 30% are needed in the ISRE intensity (up to ~ 
60% for the full range of observations). 

Acknowledgements. FB and GA acknowledge support from the ERC grant no. 
267934. 
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Fig.A.l. Scatter plot of A^/^v and fitted dust temperature for 
the modified black-body fit. Two alternative estimates of Ay are 
shown {see text), plus an average value. The horizontal grey lines 
show A^/Ay = 1 and the 15% uncertainty in the Ay normalisa¬ 
tion (Sect. 2.1). All points have similar relative uncertainty; the 
bars for the lowest are shown. The uncertainties are sys¬ 
tematic, so they affect each point in the same way; errors on the 
two axes are strongly anti-correlated. The Umin colour scheme 
is the same as Fig. 2; the filled symbols correspond to the SEDs 
shown in fig. 1 . 
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Appendix A: Modified biack-body modeis 

Fig. A. 1 shows two different estimates of Ay as a function of 
temperature, obtained from a modified black-body fit to the 
observed SEDs. Planck Collaboration XI (2014) provides two 
different empirical relations for computing the dust reddening 
E{B -V) from the emission. One relation considers E{B -V) to 
be proportional to the 850 fim dust optical depth {E{B - V)/to - 
1.49 ■ 10'^), the other relation considers E{B - V) proportional to 


the radiance (E(B - V)/'R - 5.4 ■ 10^). We calculate E(B - V) in 
both ways, and then we convert the values to Ay using the aver¬ 
age diffuse ISM value for/?y = E{B - V)IAy = 3.1 and compare 
the results. 

The two estimates differ both by their average value and their 
trend with temperature; the Ay obtained from tq is about 20% 
too low on average and decreases with temperature like the phys¬ 
ical models (Fig. 2); the Ay obtained from 'R has a good aver¬ 
age value but increases with temperature. By way of comparison 
we also show the geometric average of the two. This average 
matches the expected value better, despite having no physical 
justification. Interestingly, the two estimates implicitly make op¬ 
posite assumptions: the Ay obtained from tq assumes that the 
dust optical properties are fixed, or at least ro/Ay is fixed; and 
the Ay from R assumes a fixed Go, or at least a fixed absorbed 
power per grain . This means that cold dust is more emis¬ 
sive than expected from models with fixed dust properties, but 
less emissive than expected from models where variable optical 
properties account for all observed variations. 

Appendix B: Effects of grain size distribution 

The size distribution of dust grains has a strong effect on extinc¬ 
tion, but it is not expected to affect the far-infrared opacity, which 
only depends on the total volume of the grains in the Rayleigh 
regime. It is natural, therefore, to consider variations in grain 
size distribution as a way of varying t fir {Ay. Assessing the ef¬ 
fect of grain size variation in a physically realistic way is not 
straightforward: the physical processes that change grain sizes, 
e.g. shattering, sputtering, accretion, coagulation, also affect its 
structure and composition; also, the optical properties of the ma¬ 
terials themselves may be size-dependent. Simply varying the 
grain size distribution in a model is therefore not likely to mimic 
the actual variations in the ISM, but can still provide interest¬ 
ing qualitative insights. In this Appendix we explore the effects 
of varying the size distribution using the Cl 1 model. While this 
model does not fit the average hjAy as well as J13, it is still 
close enough to be useful for a differential analysis. Its homoge¬ 
neous grains and constant optical properties allow us to modify 
the grain size distribution independent of optical properties. 

In Cl 1, grains larger than ~ 10 nm are distributed according 
to a power law - n{a) oc where a is the grain radius - with an 
exponential cutoff above ~ 150 nm. The parameter that mainly 
controls the size distribution is the exponent of the power law, a, 
which is -2.8 for carbonaceous grains and -3.4 for silicate grains. 
We repeat the procedure of Sect. 3.2 varying a by -0.5 and -1-0.5 
around its standard value. This changes Ry by -0.7 and -i-l.O, 
respectively; by comparison, (Fitzpatrick & Massa 2007) give 
~ 0.3 as the typical 1-cr dispersion of Ry in the diffuse ISM. 

The results are shown in Fig. B.l. Varying the size distribu¬ 
tion has a small impact on the dust SED, and the range of temper¬ 
atures reproduced is smaller than that observed despite the large, 
possibly overestimated, span in a. The figure also shows that 
models with smaller grains are, surprisingly, colder than mod¬ 
els with larger grains, i.e. they have lower lOO-jum emission and 
higher long-wavelength emission. This can be explained simply. 
In our modelling, the radiation field intensity is not fixed, but 
is rather derived from the observed radiance per unit extinction 
RjAy. Changing the size distribution modifies our estimate for 

The absorbed power depends not only on Go, but also on dust prop¬ 
erties such as albedo and size distribution. Our conclusions about vari¬ 
ations in To, however, do not change if in the following we consider Go 
to mean an “effective Go” that includes albedo variations. 
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Fig.B.l. Effect of the grain size distribution on 1,1/% as pre¬ 
dicted by the Cl 1 model. The SEDs shown are the coldest (blue 
dot-dashed line), the median (green solid line) and the warmest 
(red dotted line); normalisation is the same as per Eig. 4. The cor¬ 
responding observations are plotted in grey behind the models. 
Larger symbols indicate larger average grain size, and smaller 
symbols indicate smaller average grain size (see text for details). 


Go, so that models with smaller grains necessitate a weaker ra¬ 
diation field to satisfy those constraints. The decrease in Go thus 
offsets the temperature increase due to size effects, and even re¬ 
verses it in the case of the Cl 1 model. 

The details of the result presented in this appendix are likely 
to depend on the dust model and parametrization used. Still, 
varying the grain size distribution without its corresponding 
change in the dust optical properties is not likely to explain the 
observed variations of /i/Ay. In Ysard et al. (2015) a similar 
study uses the J13 model, which is instead adapted to reproduce 
the interplay of grain size and optical properties. 
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